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OO ' ABSTRACT 

-xt : 

\l ' Aims. The kinematics, shaping, density distribution, expansion distance, and ionized mass of the nebula Hen 2-104, and the nature of its 
| symbiotic Mira are investigated. 

£f} . Methods. A combination of multi-epoch HST images and VLT integral field high-resolution spectroscopy is used to study the nebular 
' dynamics both along the line of sight and in the plane of the sky. These observations allow us to construct a 3-D spatio-kinematical model 
OO of the nebula, which together with the measurement of its apparent expansion in the plane of the sky over a period of 4 years, provides the 

■ expansion parallax for the nebula. The integral field data featuring the [Sn] AA67 1.7,673.1 emission line doublet provide us with a density map 
' of the inner lobes of the nebula, which together with the distance estimation allow us to estimate its ionized mass. 

• i-H . Results. We find densities ranging from n f =500 to 1000 cnT 3 in the inner lobes and from 300 to 500cm -3 in the outer lobes. We determine an 
' expansion parallax distance of 3.3±0.9 kpc to Hen 2-104, implying an unexpectedly large ionized mass for the nebula of the order of one tenth 
' of a solar mass. 
_ _ i 
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1. Introduction 

Symbiotic stars are valuable laboratories for studying the 
physics of a variety of astrophysical processes, like mass trans- 
fer, accretion and outflow, thermonuclear runaways, photo- 
ionization and shock-ionization, and the physico-chemistry of 
gas and dust in complex environments containing both hot and 
cold stars. 
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The formation of the nebulae that surrounds some of them, 
especially those containing a Mira variable, is one example 
of these complex interactions, and the key to understand cru- 
cial aspects of the mass loss history and geometry from these 
systems. A small fraction of the Mira wind is believed to be 
accreted by the white dwarf companion (Kenyon & Webbink 
|1984l l, resulting in quasi-stable thermonuclear burning at its 
surface (e.g. Sokoloski 2002) which sometimes ends in ther- 
monuclear outbursts (e.g. Kenyon 1986, Livio et al. 1989) 
lasting for hundreds of years (e.g. Munari 119971 ). The fast 
winds from the white dwarf collide with the slower, non- 
accreted Mira wind, and shape complex circum-binary nebu- 
lae. Currently, about a dozen of these nebulae have been opti- 
cally resolved (Corradi 2003 ). Most of them show bipolar mor- 
phologies, similar to those of many planetary nebulae (PNe) 
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and protoplanetary nebulae (PPNe) (Corradi & Schwarz 1995 
Sahai 120021 . 

The distance of these system is generally poorly known, in 
spite of being a basic parameter to determine their overall phys- 
ical properties, like their total luminosities, the mass of their 
ejecta or their kinematical ages. Recently, the advent of high 
spatial resolution instruments like the HST allows to face this 
problem by measuring the expansion parallax of some of these 
nebulae. The procedure is to measure the growth of the neb- 
ula in the plane of the sky via multi-epoch imaging, and com- 
bining this information with the expansion velocity, once the 
Doppler shifts have been measured and deprojected according 
to the adopted 3-D morphology. It has been applied to several 
relatively close PNe (NGC 6543, Reed et al. [T9991 IC 2448, 
NGC 6578, NGC 6884, NGC 3242 and NGC 6891, Palen et 
al.|2000) and to one symbiotic nebula (Hen 2-147, Santander- 
Garcfa et al. l2007al) . 

Named the Southern Crab after its resemblance to the legs 
and claws of that animal, the nebula around Hen 2-104 (also, 
PN G3 15.4+09.4, and V852 Cen) was discovered by Schwarz 
et al. (|1989l l. Although similar, in morphology and kinemat- 
ics, to several bipolar planetary nebulae (Corradi & Schwarz 
I1993bl) . the system has been unequivocally recognized as a 
genuine symbiotic star after Whitelock d 19871 ) detected a Mira 
in its core by means of long-term photometric monitoring. The 
dynamics of this magnificent nebula, which spans 70"x35" in 
the plane of the sky, has been studied by Corradi & Schwarz 
(1993a) and then analysed in great spatial and spectral detail 
by Corradi et al. J200TT ). 

The distance to Hen 2-104 is very controversial. While 
Wright and Allen ( 1978 ) found a distance of 7.92 kpc, Schwarz 
et al. (fT989l estimated 0.8 kpc. Whitelock JT987I ) and then 
Corradi et al. (2001 ) used the Period-Luminosity relation for 
Mira stars to derive an intermediate distance of 4.7 kpc and 
4.4 kpc, respectively. 

In this paper, we present the results of the analysis of an 
NTT+SOFI NIR spectrum, integral-field high-resolution VLT 
spectra, as well as of an 2003 HST [Nn] image of the nebula 
which is compared to the one obtained on 1999. With these, 
we present the first spectroscopic detection of the Mira, the de- 
termination of the expansion parallax of the nebula, its density 
map, and an estimation of its ionized mass and total luminosity 
of the system. 

2. Observations and Data reduction 

2.1. Imaging 

Two sets of [Nn] 658.3 nm images of Hen 2-104 were ob- 
tained using the Wide Field Planetary Camera 2 (WFPC2) on 
the Hubble Space Telescope (HST) on May 22, 1999 (pro- 
gramme G07378) and May 10, 2003 (programme G09336) 
respectively (see Fig. Q~|(. Images from both epochs were ob- 
tained through the F658N filter (central wavelength/bandpass = 
659.0/2.2 nm), which, assuming no significant degradation be- 
tween epochs, transmits the [Nn] line with 78% efficiency and 
the Hot line with 1 % efficiency at the heliocentric systemic ve- 
locity of Hen 2-104. An additional set of WFPC2 images was 



taken in May 20, 2001 (programme GO9050), in the F656N fil- 
ter (central wavelength/bandpass = 656.4/2.2 nm). In this case, 
the transmission efficiency is 78% in Ha, 11% in [Nn] 654.8 
and 6% in [Nn] 658.3. 

We computed the observed total fluxes from the nebula in 
both Ha and [Nn] 658.3 nm following the procedure described 
by Luridiana et al. (2003 ). After subtraction of the central un- 
resolved core and decontamination from continuum and neigh- 
bour emission lines, the observed total fluxes from the nebula 
were estimated to be 3.6 10~ 12 in Ha and 2.1 10~ 12 erg ctrT 2 s _1 
in [Nn] 658.3 nm. These are consistent with the values given by 
Shaw & Kaler (119891 ) for a 10" diameter aperture centered on 
the star. 

The HST images (one dataset per epoch) were obtained us- 
ing a 4-point dithering pattern in order to properly sample the 
HST point spread function (PSF). They were then reduced us- 
ing the "multidrizzle package" (Koekemoer et al. 2002}, an 
automated dithered image combination and cleaning task in 
PyRAF (Python IRAF). This routine first creates a bad-pixel 
mask using calibration reference files, and performs the sky 
subtraction. Next, the software arranges the four WFPC2 CCDs 
into a mosaic, where geometric distortion effects are corrected, 
including time dependent distortions due to CCD rotations and 
flexure. Then it removes cosmic rays and creates, for each 
epoch, a combined and drizzled image with a sampling scale 
of0".045pix- 1 . 

In order to extract the expansion on the sky of Hen 2-104 
from these two images, they need to be brought to exactly the 
same astrometric grid. This is not the case due to slightly dif- 
ferent roll angles and pointings, and possible residuals in the 
distortion correction done with multidrizzle for the two epochs. 
We used THEL.£] (Erben et al. 1995 ) to define an astrometric 
reference frame based on the stellar sources in the field. Given 
the small rotation and distortion corrections, the images were 
then resampled to the new grid, using the coordinates of the 
central star of Hen 2-104 as the reference projection point. The 
relative astrometric accuracy between the two images is on the 
order of 0.2 pixel. 

2.2. Integral Field Spectroscopy 

Integral Field Spectroscopy was obtained at Kueyen, Telescope 
Unit 2 of the 8 m VLT at ESO's observatory on Paranal, on 
April 2, 2004. The inner lobes of the nebula were covered by 
the ARGUS IFU at PA=122° using the lowest spatial resolu- 
tion (314 fibres each with a size of 0".52 projected onto the 
sky, giving a total of 11". 4 by 7". 3 field of view). Two 300 s 
exposures were taken with the H679.7 grating, covering the re- 
gion between 660 and 690 nm, and then co-added. The resolv- 
ing power was R=30000, corresponding to a line FWHM of 
10.5 km s _1 . The seeing was 0".7. The VLT data were kindly 
reduced by Reinhard Hanuschik through the ESO pipeline, and 
then further analysed using standard IRAF routines. 
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Fig. 1. The HST F658N image of the nebula of Hen 2-104 obtained 
in 2003. The ARGUS IFU unit covering the inner lobes is overlaid 
on the top-left zoomed region. 
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Fig. 2. NTT+SOFI NIR spectrum of the core of Hen 2-104. Regions were telluric absorption bands severely degraded the 
spectrum have been removed. Right bottom: Magnification of the region spanning from 2.2 to 2.45 jjm. 



2.3. Near-infrared spectroscopy 

A 120 s near-IR spectrum of the core of Hen 2-104 was taken 
with SOFI at the ESO New Technology Telescope (NTT) on 
May 2, 2004, in the 0.95-2.5 yum range, with a resolving power 
of R-930 for a slit width of 0".6. The spectrum was reduced 
following standard IRAF procedures. 

3. An improved spatio-kinematical modelling 

The image of the nebula surrounding Hen 2-104 is shown in 
Fig. [1] It consists of two pair of nested, emitting gas lobes. 
According to Corradi et al. (2001), the two pairs of lobes are 



likely coeval, and not the result of two similar consecutive out- 
bursts (with the inner lobes being a younger version of the outer 
ones) since their expansion velocities are (dramatically) differ- 
ent, and their kinematical ages are very similar. Also the polar 
jets of the nebula are roughly coeval with the inner and outer 
lobes. The inner lobes show latitude dependent surface bright- 
ness variations, which produce the systems of ellipses observed 
in the HST images (Corradi et al. 2001). The outer lobes, on the 
other hand, appear fragmented by Kelvin-Helmholtz instabili- 
ties into a myriad of clumpy and filamentary structures. Most 
of them, specially in the furthest regions, have outward fac- 
ing tails (see Contini & Formiggini 120011 Gonzalez, Raga & 
Steffen|2j05}. 
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We used the VLT IFU to map the 2-D distribution and 
Doppler shifts of the emission from the core and the whole in- 
ner lobes. The core shows a bright Hei /1667.8 emission line, 
which is also visible throughout the inner nebula. The [Sn] 
/Lt671.7,673.1 emission lines present a high S/N ratio over the 
whole inner lobes, while are somewhat weaker in the portion 
of the outer lobes contained in the small field of view of the 
IFU. This, together with a wavelength dispersion large enough 
to separate the emission from the two system of lobes wherever 
their respective radial velocities differ in more than 40 km s , 
allows us to make a 3-D analysis of the two structures sepa- 
rately. 

The IFU [Sn] 22x14 fiber-array was rearranged in columns 
and rows to construct a grid of pseudo long-slit spectra, with 
length 11. "4 and 7. "28 respectively, and slit width 0."52. A 
spatio-kinematical model was then built and fit to the spectral 
data and to the shape of the inner lobes and base of the outer 
lobes in the [Nn] HST image (see Santander-Garcfa et al. 2004 
for a detailed description of the method). An analytical descrip- 
tion assuming axial symmetry and using the formula in Solf & 
Ulrich ( 1985 ) was adopted, with a surface, described in spher- 
ical coordinates by: 



equator 



+ (V polar ~ Vequator) SU1 |6f] 



where r is the apparent distance to the central source, tD' 1 
the kinematical age of the outflow per unit of distance to the 
nebula, v po iar an d ^equator the velocities of the model at the pole 
and equator respectively, 6 the latitude angle of the model, and 
y a dimensionless shaping factor. This assumes that each gas 
particle travels in the radial direction from the central star with 
a velocity proportional to its distance to the central source. For 
the inner lobes, emission is also restricted to the latitudes that 
correspond to the observed system of ellipses. In both cases, 
the resulting two-dimensional model was scaled to fit the size 
of the object, rotated around the symmetry axis to produce a 
three dimensional figure, and inclined to the plane of the sky. 
The resulting geometrical shape and velocity field is used to 
generate simplified model images and velocity plots for direct 
comparison with the real image and spectra. 



Parameter 


Value 


Range 


tD [ (yrs kpcr 1 ) 


1300 


(1200-1400) 


VpoUr (km S _1 ) 


230 


(210-250) 


Vequator (km S ) 


12 


(8-15) 


y 


3.2 


(2.8-3.8) 


Omax (°) 


52 


(50-53) 


in 


58 


(56 - 60) 



Table 1. Best-fitting parameters for the outer lobes of Hen 2- 
104. 



The fit to the data is performed by allowing the different pa- 
rameters of the model to vary over a large range of values and 
visually comparing each resulting model to the image and spec- 
tra, until the optimal model is found. The results are shown in 
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Fig. 3. Left: The [Nn] HST image of the inner and outer lobes 
of Hen 2-104. Right: Adopted model of the outer lobes. 




Fig. 5. Left: The [Nn] HST image of the inner lobes of Hen 2- 
104. Right: Adopted model of the inner lobes. 



Parameter 


Value 


Range 


tD [ (yrs kpc" 1 ) 


1100 


(900-1300) 


v polar (km s _1 ) 


90 


(210-250) 


Vequator (km S ) 


11.3 


(7-15) 


r 


4.5 


(3.5-5.0) 


@min ( ) 


25 


(24 - 27) 


$max ( ) 


39 


(37 - 40) 


i(°) 


58 


(56 - 62) 



Table 2. Best-fitting parameters for the inner lobes of Hen 2- 
104. 



Figs. [3] |U [5]and[6]and the corresponding parameters in Tables 
1 and 2. 

Almost all the emission from the inner lobes, both in the 
[Nn] image and the [Sn] spectra, was found to be faithfully 
accounted for by a series of rings inscribed on the walls of a 
truncated Solf & Ulrich ( 119851 model at the highest latitudes, 
from 25° to 39°, and expanding at a maximum velocity of ~25 
km s . The outer lobes, on the other hand, were found to share 
the symmetry axis of their smaller counterpart, 58±2° inclined 
with respect to the line of sight, and with a similar, within un- 
certainties, age of the nebula. We thus confirm the overall re- 
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Fig. 4. The ESO [Sn] ARGUS IFU spectra rearranged so as to simulate long-slit spectra at P.A. = 122° (long axis of the nebula), 
and offset from -3". 38 (SW) to +3". 38 (NE) from left to right, in increments of 0".52. The spatio-kinematical model of the outer 
lobes has been superimposed. Each frame is 10". 92 tall and 260 km s _! wide. Southeastern side is up. 



( 'i i t i m * "i 



Fig. 6. The ESO [Sn] ARGUS IFU spectra arranged in columns along the IFU longer axis. Frames correspond, from left to right, 
to P.A.=122° and offsets from -3". 38 (SW) to +3". 38 (NE) from left to right, in increments of 0".52. The spatio-kinematical 
model of the inner lobes has been superimposed. Note: Each frame is 10". 92 tall and 260 km s _1 wide. Southeastern side is up. 



suits by Corradi et al. (2001), while refining the value of the 
fitting parameters for the inner lobes thanks to the complete 
coverage of this region of the nebula by the ARGUS IFU. A 
systemic velocity \lsr = -85 km s was derived, in agreement 
with the measurement of Corradi & Schwarz dl993al >. 

4. The distance 

The distance to the nebula of Hen 2-104 has been determined 
by means of its expansion parallax. This method calculates the 
distance by computing the apparent expansion in the plane of 
the sky between two images taken in two different epochs, and 
combining it with a 3-D model of the expansion pattern, in turn 
usually recovered from the analysis of the images and the radial 
velocity field. 

The apparent expansion of both the outer lobes and po- 
lar jets of Hen 2-104 is immediately visible by blinking the 
aligned 1999 and 2003 HST images. The amount of expan- 
sion of the outer lobes and jets was quantified in two different 
ways (Reed et al. 1999): 1) deriving the growth of the whole 
outer lobes and jets with the so-called magnification method; 
and 2) working out the expansion of several individual knots 
via cross-correlation of the corresponding surface brightness 
radial profiles. 

4.1. Expansion parallax via the magnification method 

A residual image of the nebula was obtained by subtracting 
the 1999 image from the 2003 one, once the images have been 



registered and aligned with respect to the central core. A gen- 
eralized pattern can be seen in the jets and edges of the outer 
lobes, with positive residuals on the outside and negative ones 
on the inside, indicating the expansion of the nebula in the pe- 
riod from 1999 to 2003. 

In order to quantify the growth of the nebula, we deter- 
mined the magnification factor, M, to be applied to the 1999 
image, which minimizes the rms of the residual image ob- 
tained by subtracting the magnified image from the 2003 im- 
age. In practice, though, the complex structure of the nebula, 
the evolution of the brightness and shape of individual knots 
and the small, non-radial components of movement of each 
knot (Rayleigh-Taylor instabilities, diffusion into the medium, 
etc.) complicate the analysis. 

Each jet, as well as two well exposed regions of the outer 
lobes (see Fig|7]i, were independently analysed. The minimum 
rms of the residual images of the Northwest and Southeast jets 
were found to correspond to magnification factors M equal to 
1.0009 and 1.0010 respectively; both with an uncertainty of 
+0.0002. On the other hand, the M factors which minimized 
the residual images of the Northwest and Southeast regions of 
the outer lobes were 1.0008 and 1.0012 respectively, while the 
uncertainty in these cases was +0.0003. Independent compar- 
isons of the residual images by eye confirmed all these results. 
Thus, the growth of both the jets and outer lobes is uniform, 
within the errors, with average values of M-l. 00095+0.00020 
and M=1.00100±0.00030, respectively. 
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Fig. 7. HST [Nn] WFPC2 2003 image. The four regions in which 
the rms of the magnification method residual image were measured 
are indicated by the black boxes. 



Given the magnification factor M for the jets and outer 
lobes, and the kinematical age-distance parameter tD 1 (1300 
yrs kpc~' for the outer lobes, from Tab. [3] and 1200 yrs kpc 1 
for the jets, Corradi et al. |T993a), the expansion parallax can 
be determined by comparing the amount of growth of each fea- 
ture in the lapse of time between 1999 and 2003, tg, with the 
features' kinematical age. Hence: 



D(kpc) = 



tsiyrs) 



(M - 1) tD- l {yrs kpc- 1 ) 



We then found that the magnification method applied to the 
jets and outer lobes of Hen 2-104 gives consistent values of 
Dj ets -3. 5+0. 7 kpc and D outer i t, es -'i. 1+0.9 kpc, respectively. 



the cross-correlation of the radial profile. The inclination is as- 
sumed to be the same as for the inner and outer lobes (Sect. 3). 
For the outer lobes, the velocity in the plane of the sky is com- 
puted according to the spatio-kinematical model in Sect. 3. 

After averaging the results for all the knots in the outer 
lobes, we obtained D OMer i \, es = 3.1 + 1.0 kpc, whereas Dj ets = 
3.4+ 1 .0 kpc for the knots in the jets. 

The distance found through radial fitting is therefore 
fully consistent with the one found through the magnification 
method; in the following, we will adopt an average distance 
of Duen 2-104 = 3.3+0.9 kpc. Note the considerable size of the 
nebula at this distance, amounting to ~2.0 pc, and the large 
kinematical ages of the inner and outer lobes, which would be 
3700+1200, and 4200 +1200 years, respectively. 



4.2. Expansion parallax via radial profile method 

While not providing a global, intuitive vision of the growth of 
the nebula and of deviations from pure radial expansion, the 
method that involves the study of radial profiles allows us to 
trace the radial expansion of individual features. 

Essentially, this procedure consists of the analysis of one- 
dimensional surface brightness radial profiles passing through 
a given bright and sharp knot. Ideally, any profile from the 2003 
image will be similar to the one from 1999, except that the fea- 
ture will be shifted outwards. In order to measure the angular 
expansion 6 of the feature in the plane of the sky, the 1999 and 
2003 profiles are cross-correlated. 

Fig. [8] shows the angular expansion 6 (in milliarcseconds) 
of 45 well-defined knots chosen in the outer lobes and jets. 

Assuming pure radial motion of the features, the distance 
to the nebula is computed by comparing the angular expansion 
of each knot to its velocity in the plane of the sky: 



D(kpc) = 0.211 



v sky (km s ') 



For each knot in the jets, its radial velocity from the ma- 
jor axis spectrum at P.A.=122° from Corradi et al. (2001), was 
coupled with the expansion in the plane of the sky derived from 



5. The NIR spectrum 

The NIR spectrum of the core (see Fig. [2]) shows a variety of 
emission lines over an almost featureless rising continuum. 

The region between 0.95 and 1.3 //m shows Hydrogen and 
Helium atomic recombination lines, as well as [SII], [SHI], 
[NI], and H2. The region onward 2.1 yum shows three CO ab- 
sorption bands typical of late AGB stars (e.g. Kamath & Ashok 
2003). These bands, although barely visible in the dominant 
continuum presumably produced by hot circumstellar dust, 
show for the first time some photospheric spectral features of 
the Mira in the system, so far only detected thanks to its char- 
acteristic NIR photometric light curve (Whitelock l 19871 ). 

The emission lines and absorption band CO identifications 
are shown in Table 3. 

The featureless NIR continuum of Hen 2-104 is quite pe- 
culiar for symbiotic stars (cf. Whitelock 2003 ), where in gen- 
eral this part of the spectrum is dominated by deep absorption 
bands of the red giant. Note that Mz 3 and M 2-9, two nebulae 
suspected of having a symbiotic binary nucleus, show a simi- 
lar spectrum. It is likely that the rising continuum is emission 
produced by a significant amount of hot circumstellar dust, but 
a detailed modelling of the spectrum is beyond the scope of 
these paper. What is instead relevant for the present discussion, 
is the fact that it is clear that the emission in the K-band is not 
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Fig. 8. HST [Nil] WFPC2 2003 image along with the radial expansion (in mas yr 1 ) of the knots of the outer lobes and jets 
analyzed in section 3. 



always a good measure of the luminosity of the red giant of 
the system, so that distances computed from the Mira period- 
luminosity relationship using the K magnitude, are necessarily 
subject to large uncertainties. 



cm -3 is estimated for Hen 2-104 from the [Om] lines (Freitas 
Pachecho & Costa 1996). At these densities, the [Sn] lines are 
collisionally quenched. These results confirm those given by 
Corradi & Schwarz (|1993al >. 



6. Nebular density and mass 

6.1. Density map 

The integral field spectroscopy data, featuring the [Sn] 
/t/1671.7,673.1 emission line doublet, allow to build a den- 
sity map covering the whole inner lobes and the base of the 
outer lobes (Fig |9j . As mentioned previously, the large resolv- 
ing power (R=30000) of the spectra was high enough to sepa- 
rate the inner lobes from the outer lobes wherever their respec- 
tive radial velocities differed more than 40 km s _1 . 

In the absence of actual temperature measurements for the 
lobes, we assume T e - 10000 K for both structures. The result- 
ing densities range from n e =500 to 1000 cm -3 for the inner 
lobes, and from ~300 to 500 cm -3 for the base of the outer 
lobes. The core shows the highest [Sn] densities, n e ~2400 
cm -3 . However, the [Sn] lines seen towards the core probably 
arise in the outer layers of a very dense circumstellar envelope, 
as e.g. in the case of the other symbiotic Mira HM Sge (Corradi 
& Schwarz Tl999l l. Indeed, a core density larger than 2.5 10 6 



6.2. Ionized mass 

The ionized mass of the nebula around Hen 2-104 can be esti- 
mated via the formula 

m H F(H/3)47TD 2 
Mh+ ~ u ~f 

" v /3 a B n <? 

where D is the distance, n e the nebular electron density, 
F(Hy6) is the total H/3 flux, niu is the mass of the hydrogen 
atom, hvp the energy of an H/3 photon, and the effective 
recombination coefficient for case B (adopted from Osterbrock 
1989 ). After correcting the observed Ha flux for an extinction 
of A,,=3.0 mag (average of the values found by Lutz et al. 1989 
and Freitas Pacheco & Costa 1996), and applying the theo- 
retical Ha to Hy6 ratio of 2.85 (for an electron temperature of 
10000 K), we obtain F(H/3)=1.3 lO" 11 erg cm 2 s -1 . Assuming 
an average density of 500 cirT 3 and the distance found in this 
work, we compute an ionized mass of ~0.16 M Q . The extinc- 
tion correction is very important in the above calculation, and, 
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A Q (laboratory) Element Flux error 

(A) (emission) (1(T 13 (%) 

14 
1 

17 
1 
1 

10 

7 

20 
7 
1 

18 
1 

10 
31 
36 
18 
10 
10 
7 
1 

39 
52 
15 



A (bandhead) 


Molecule 


Flux 


error 


(/an) 


(absorption) 


(10 -13 «!_■) 


(%) 


2.32 


CO 


0.34 


28 


2.37 


CO 


0.35 


46 


2.42 


CO 


0.59 


27 



Table 3. Hen 2-104: NIR emission line and absorption band 
identification, along with their fluxes and associated errors, 
which include both the poissonian error and the error in the 
determination of the continuum. "?" denotes that the line was 
unidentified. 



if we instead assume an interstellar extinction of 1 .6 mag, as de- 
rived from the Galactic extinction model by Drimmel, Cabrera- 
Lavers & Lopez-Corredoira (2003) for the location of Hen 2- 
104 and the computed distance, the resulting ionized mass is 
0.05 M . These have to be considered as lower limits, as they 
neglect any local extinction due to dust around and inside the 
nebula. 

A different mass estimate can be obtained by computing, 
both for the inner and outer lobes, the volume of each feature 
occupied by gas in the following manner. As the geometry of 
the surface of the lobes is known from our spatio-kinematical 
modelling (Sect. 3), and the electron density distribution is 
also known from the [Sn] line ratio, the only other parame- 
ter needed to compute the volume occupied by the emitting gas 
is the thickness of the walls of the lobes, if the filling factor 
(cf. Pottasch 1984) is assumed to be around unity. 3-D models 
with different thickness have been produced using IDL (RSI) 
routines adapted from the VISNEB tools developed for the 
Cloudy 3D code (Morisset l2006b . A cubic grid of 300 3 cells 



J L 




Fig. 9. Density map of the inner lobes and base of the outer 
lobes of Hen 2-104. The contour levels span from n e - 200 to 
3000 cm 3 with step 200. 

was used to generate and rotate a 3D emissivity model (with 
the emissivity per unit volume oc n e 2 ). The emission along the 
line of sight was then integrated in order to generate a synthetic 
image in collisionally excited lines (this also assumes that the 
ionic abundance and electron temperature are constant through- 
out the lobes). This procedure was repeated several times al- 
lowing the thickness of the 3D model to vary over a wide range. 

In order to find the optimal thickness that accounts for 
the observed inner and outer lobes, residual images were ob- 
tained by subtracting each synthetic image from the real HST 
one. The regions around the edges of each lobe (i.e. where the 
walls of the structure are perpendicular to the plane of the sky 
and therefore its thickness is easier to estimate) were anal- 
ysed in the residual images, both by visual comparison and 
by measuring the rms of each region. Illustrative samples of 
these synthetic images, both for the outer and inner lobes, are 
shown in FigsJlO] and QT| The resulting best-fit thickness is 
0".9+0".4 and 3".4+0".5 for the inner and outer lobes respec- 
tively. Finally, in order to account for the knotty appearance of 
the outer lobes, we computed the fraction of the total volume 
of the model occupied by gas. The procedure to do this was 
to count the number of pixels in the HST images whose sur- 
face brightness is larger than 3 times the statistical noise of the 
background, and compute which fraction of the area projected 
in the sky of the model nebula is effectively filled. The result 
was that 0.34+0.06 of the total volume of the outer lobes is 
filled by gas. 

Therefore, using typical densities of n e =600 cm -3 and 
n e =350 cirT 3 for the inner and outer lobes, respectively, and 
assuming a roughly solar abundance for the nebula, the result- 
ing ionized mass for the two pairs of lobes is (4.6+0.9) 10~ 3 
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and 0.08+0.01 M , respectively, for a total of roughly ~0.1 M 
for the whole nebula. This value is midway, but broadly con- 
sistent, with the two previous estimates, 0.16 and 0.05 M Q ob- 
tained varying the adopted extinction, and is clearly pointing to 
a relatively large mass for the nebula. 

7. Discussion 

Schwarz et al (1989) discovered the extended nebula around 
Hen 2-104, and they would probably have classified it as a 
planetary nebula, had Whitelock (1987 ) not found, two years 
before, a long-period pulsation of the central star by means of 
NIR photometric monitoring. In fact, apart from indirect hints 
suggesting a symbiotic nature (e.g. a high density core, NIR 
colours unusual for a PN, and an articulated high-excitation 
spectrum), this was the only evidence for the presence of a cool 
giant. The NIR spectrum presented in Sect. 5 shows for the 
first time the spectral signature of the red giant in the spectrum. 
However, a strong featureless continuum, presumably from hot 
circumstellar dust, renders these bands barely visible. Thus, 
given the difficulties of direct detection of the Mira in some 
genuine symbiotic nebulae such as Hen 2-104, the already long 
debate on the nature of some bipolar nebulae with dust domi- 
nated NIR spectra and high density compact cores (e.g. Mz 3 
and M2-9, see Schmeja & Kimeswenger I20031 Phillips [20071 
and Santander-Garcfa & Corradi 2007b) is not a surprise. 

In this work, we derive a distance for Hen 2-104 of 
3.3+0.9 kpc by means of the expansion parallax method. It is 
known that the expansion parallax determinations can be af- 
fected by shocks. This is for example the case of the symbi- 
otic Mira Hen 2-147, in which Santander-Garcfa et al. (2007a) 
found that distance is understimated if shocks are not properly 
taken into account. The main evidence for shocks in that neb- 
ula are the broad and asymmetrical line profiles, amounting to 
up 200 km s _1 full width at zero intensity (cf. Hartigan et al. 
119871 ). This is not the case for Hen 2-104, where line widths in 
the jets and lobes are of the order of ~15 km s _1 , only slightly 
larger than the gas thermal broadening. We therefore conclude 
that shocks do not strongly affect the expansion parallax deter- 
mination of Hen 2-104, as also suggested by the consistent dis- 
tance estimate for the lobes and the jets (shocks would a priori 
act differently on the two components Jl We consider our ex- 
pansion parallax determination to be more robust than the dis- 
tance obtained through the period-luminosity relationship for 
the Mira (Whitelock 1 1 9871 ), because our spectra show that the 
K band magnitude does not provide a reliable measure of the 
luminosity of the Mira. 

Application of period-luminosity relation (Feast et al. 1989 
Groenewegen & Whitelock [l9961 and Whitelock et al l 19941) to 
the Mira of Hen 2-104 (P=400 days, Whitelock [19871 , would 
imply a bolometric luminosity of the AGB star of ~6000 L Q . 
This can be compared with the total luminosity of the sys- 
tem derived from the spectral energy distribution presented 
in (Fig. IT2b . at the computed distance of 3.3 kpc. From the 

2 Note, however, that even if the actual distance were somewhat 
larger, the conclusions below would still hold, as the ionized mass 
would increase accordingly. 
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Fig. 12. Observed spectral energy distribution of Hen 2-104, 
from 1 .2 A to 100 fim, based on the values quoted in Schwarz et 
al. [19891 for the UBV and the IR photometry, the IRAS fluxes, 
the IUE UV fluxes and the Chandra X ray fluxes. 



latter, in order to obtain a luminosity larger than 6000 L Q 
(the luminosity of the white dwarf companion should also be 
considered), the interstellar extinction must be A v >2.9 mag. 
This value is significantly larger than the one obtained us- 
ing the Galactic model by Drimmel, Cabrera-Lavers & Lopez- 
Corredoira ( 120031 see Sect. 6.2), and closer to the determina- 
tions by Lutz et al. (1989), A,,=2.23 mag, and Freitas Pacheco 
& Costa ( [19961 ), A v ,=3.7 mag. 

Another result of this work that deserves some discussion is 
the remarkable ionized mass of the extended nebula, amount- 
ing to about a tenth of a solar mass (see section 6.2). Given 
such a large value, and considering that only a negligible frac- 
tion of the Mira wind is captured by the white dwarf, it seems 
clear that the main donor of mass for the nebula is the Mira. 
As Corradi et al. (2001) suggested, this kind of nebulae are 
likely to be the result of the interaction of the slow winds from 
the Mira with the fast winds from an outbursting white dwarf. 
Given the articulated structure of these nebulae, it is not easy to 
understand how the different winds contribute to each morpho- 
logical component, but it seems natural to associate the high ve- 
locity features (like the polar jets) directly with the fast winds 
from the hot white dwarf and its accretion disc. Many of the 
features of the nebula around Hen 2-104 could be explained 
by the weak-jet model of the hydro-dynamic calculations by 
Garcfa-Arredondo & Frank (2004). Although the initial condi- 
tions they assumed for the two case studies are different, the 
conclusions of their model are applicable. This emphasizes the 
key role of the Mira in supplying the material and the essen- 
tial role of the white dwarf companion in accelerating and, in 
all probability, deflecting the streamlines of the Mira's wind to 
shape the bipolar morphology that we see today. Our spatio- 
kinematic fit clearly shows that the outflow speeds are slowest 
at low latitudes and highest along the polar direction. This is 
the same pattern as found in most bipolar PNe (e.g. Corradi 
2004), as well as in other classes of objects like in the ho- 
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Fig. 10. From Left to Right: HST [Nn] image from 2003 and illustrative sample of synthetic images of the nebula. The thickness 
of the outer lobes of the three models are 1".0, 3". 5 and 6".0, respectively, while the inner lobes show a constant thickness of 
0".9. 




4>4 



Fig. 11. From Left to Right: HST [Nn] image from 2003 and illustrative sample of synthetic images of the inner lobes (and the 
central region of the outer ones). The thickness of the inner lobes of the three models are 0".4, 0".9 and 1".4, respectively, while 
the outer lobes show a constant thickness of 3".5. 



mologously expanding bipolar homunculus of r\ Carina (Smith 
& Gehrz 1998 ). Thus, the mechanisms that may have formed 
Hen 2-104 are not limited to symbiotic stars and PNe, and may 
be common when the mass loss rate of the source is high. 

We have also found that the walls of the two lobe pairs 
of Hen 2-104 are very thin compared to their radial distance 
from the star. A comparison to the sound crossing time allows 
for some speculation about the mechanism holding the walls 
tightly together. If the gas temperature is (and has been) 10 4 K 
since the ejection ~4000 yrs ago (at the computed distance), the 
sound speed is 10 km s _1 , so it travels 1.2 10 12 km in 4000 yrs. 
The wall thickness of the inner lobes (0".9), 3.5 10 11 km, is less 
than the sound travel distance, implying that something like an 
outside pressure or an embedded magnetic field is holding the 
walls of the inner lobes together, perhaps accounting for their 
relatively smooth and structured appearance and high emission 
density. The outer walls are 1.3 10 12 km thick (3". 4), which is 
roughly the sound crossing distance. This means that the walls 
of the outer lobes are expanding at their sound speed (i.e. driven 
by their own thermal pressure). 



For a typical AGB wind speed of ~15 km s _I (Habing, 
Tignon & Tielens 1994), it would take ~2.6 10 4 yrs to the gas 
lost by the Mira of Hen 2-104 to fill the volume presently occu- 
pied by the nebula (~0.4 pc in radius at the adopted distance). 
During this lapse of time, an average mass loss rate of ~3.5 
10~ 6 M Q yr~' would be needed to accumulate the total mass of 
the nebula (~0.1 M ) around the system, that later would be 
shaped by fast winds from the outbursting white dwarf com- 
panion. This mass loss rate is in good agreement with the one 
expected from a 6000 L Mira variable (Alard et al. [2001), 
making the depicted scenario a coherent one and reinforcing 
our distance and mass estimates. 

In any case, the mass of the nebula of Hen 2-104 is re- 
markably larger than previously thought for this class of ob- 
jects, and more similar to the ionized masses of genuine PNe 
(Pottasch fl984| l. Thus the Mira in Hen 2-104 would be in a 
quite evolved stage, close to the tip of the AGB, where due to 
strong mass loss a massive circumstellar envelope can be built 
up. The ionized nebula that we observe can therefore be seen as 
a sort of "premature" PN, in which the nebular shaping and ex- 
citation is done by the WD companion, which "anticipates" the 
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post-AGB evolution of the Mira. Note that the kinetic energy 
of the nebula of Hen 2-104 (~8 10 44 erg) and its momentum 
(~10 39 g cm s -1 ) are similar to those of regular PNe, and sig- 
nificantly lower than the momentum computed by Bujarrabal 
(2001 ) for a number of multi-polar young and proto-PNe. The 
energetics of the outflow of Hen 2-104 can be accounted by 
the Mira luminosity and the energy released by a slow-nova 
explosion on the WD companion. 

In the past, it has been often stated that a distinctive char- 
acteristic of the nebulae around symbiotic stars is that they are 
one or two orders of magnitude less massive than PNe. Hen 2- 
104 shows that this is not the rule, so that the line dividing PNe 
from symbiotic nebulae might actually be more blurred than 
previously thought. 
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